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Abstract: With the Thomson scattering (TS) system in KSTAR, temporal evolution of electron 
temperature (T^) is estimated using a weighted look-up table method with fast sampling (1.25 or 2.5 
GS/s) digitizers during the 2014 KSTAR campaign. Background noise level is used as a weighting 
parameter without considering the photon noise due to the absence of information on absolute pho¬ 
ton counts detected by the TS system. Estimated electron temperature during a relatively quiescent 
discharge are scattered, i.e., 15% variation on Te with respect to its mean value. We find that this 
15% variation on Te cannot be explained solely by the background noise level which leads us to 
include photon noise effects in our analysis. Using synthetic data, we have estimated the required 
photon noise level consistent with the observation and determined the dominant noise source in 
KSTAR TS system. 
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1. Introduction 

Thomson scattering (TS) systems are widely used to measure the temperature (T^) and density (n^) 
of electrons in fusion devices. KSTAR Thomson scattering system consists of separate core and 
edge collection optics collecting the Thomson scattered photons which are transmitted to polychro- 
mators via optical fibres [HI, 0]. After the band-pass filters, where Figure [T| shows an example of the 
measured filter functions, optical signals are converted to electrical signals and recorded with the 
newly installed fast sampling digitizers (NI PXIe-5160) with the sampling rate of either 1.25 or 2.5 
GS/s, in addition to the existing charge-integrating (gating) digitizer. 

We use a weighted look-up table method to estimate the time evolution of electron temperature 
while exclusively considering the uncertainty due to the background noise. Another possibly large 
source of uncertainty may originate from the photon (Poisson) noise. However, the photon noise 
could not be reflected in our analysis due to the lack of absolute photon counts. Note that the 
absolute photon counts can be estimated based on Rayleigh or rotational Raman calibration [0,0], 
but such data have not been obtained with the fast sampling digitizers during the 2014 KSTAR 
campaign. Rayleigh calibration data are obtained only with the existing charge-integrating (gating) 
digitizers. 

Estimated electron temperature shows scatters, i.e., 15% variation on Te with respect to its 
mean value. By using the synthetic data, we have found that this scattering on Te cannot be ex¬ 
plained solely by the uncertainty due to the background noise. Our main goal of this paper is to 
estimate the ‘total’ noise level of the Thomson signal consistent with the observed variation on Te 
using the synthetic data. With the levels of estimated total noise and the measured background 
noise, we deduce the photon noise level assuming that photon noise and background noise are un¬ 
correlated. We confirm that photon noise prevails over the background one for the Thomson data 
from the 2014 KSTAR campaign. 
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Figure 1. Measured filter functions (‘+’), i.e., relative transmittance coefficients of the optical band-pass 
filters as a function of the wavelength, for five separate channels of a core polychromator. Channel 1 looks 
at the shortest wavelength range, while channel 5 contains the Nd-YAG laser wavelength of 1064 nm. Also 
shown (lines) are the spectral distributions of the scattered photons O for various temperatures. 


2. Estimating electron temperature 

A polychromator contains five channels consisting of five band-pass filters shown in Figure [T| with 
five photon detectors (IR-enhanced Si APD: Hamamatsu SI 1519). Figure @ shows an example of 
the Thomson signal from each channel. Although the background noise levels are quite similar 
for all five channels, ch 1 does not show observable Thomson scattered signal. This may be due to 
electron temperature being less than 1 keV (see Figure |T|). The signal from ch 5 is mostly due to the 
straylight as attested by the fact that the signal levels before and during plasma shots are identical 
(not shown in the figure). The level of Thomson signal from ch 4 is smaller compared to signals 
obtained from chs 2 and 3, again perhaps due to not-high enough electron temperature. 

To estimate the integral of the /th channel Thomson signal we perform a Gaussian fitting 
to the Thomson signal (red line in Figure [2(b)| ) which can be written as 0] 

V^sit)^Gnenlaser{t)^^QLT{XL)QE ( 2 . 1 ) 

where G is the APD gain factor, niaser (0 number of photons per unit time as a function of 
time where / dt niaser (0 ^ ^laser is the total number of photons in a single laser pulse. is the 
differential Thomson scattering cross-section area, AQ the solid angle of the TS system and L the 
scattering length. The net transnussion coefficient of the collective optics r(A) is a function of the 
wavelength, but we assume that it is constant within the range of interests; hence r(A) = T{Xl) 
where Xl is the laser wavelength. As the /th channel filter function includes the wavelength 
variation of the quantum efficiency QE of the APD detector, the value of QE is taken at the laser 
wavelength. S{X^Te^0) is the spectral distribution of the scattered photons [0 where G is the 
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Figure 2. Measured Thomson signal from each channel for KSTAR shot #10433 at plasma time = 10 s. 
Thomson signal from ch 1 cannot be distinguished from the noise, while the signal from ch 5 is dominated 
by the straylight. Red line in (b) illustrates the Gaussian fitting with the width of w to the raw Thomson 
signal (black). Blue shaded region is used to estimate the background noise level of ch 2. 


scattering angle. Then, is simply 


— 

/\ts — 



( 2 . 2 ) 


From Eq. & we find that the terms outside the integral do not depend on the channels. 
Thus, by taking the ratio of measured signals between the two channels such as 




(2.3) 


we can construct an equation that depends on only and 6. Note that we have 6 « 90° in this 
study [0|. Defining M' {1^) = JdX </»'(A) S{X,Te, 6 = 90°), Figure t3(a)| shows M' as a function 
of Tg for all five channels. As we decide not to use signals from chs 1 and 5 due to indiscernible 
Thomson signals from the background noise or stray light, we create with chs 2, 3 and 4 which 
are shown in Figure |3(b)| . 
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Figure 3. (a) {Te) and (b) {Te) the ratio of between channels of 2, 3 and 4. 


Because we may find three different T^’s based on the measurements from Figure |3(b)| , the T^’s 
can be averaged with a weighting factor of the inverse of corresponding total uncertainty as 

L f L ’ (;>/,n>m), (2.4) 

, 7=2 ^total \m,n=2 ^total J 

where Te^ and are the estimated electron temperature using and its total uncertainty, 
respectively, can be estimated using a propagation of uncertainty assuming that the total 

uncertainty of the /th and jth channels, and are uncorrelated. Ideally, needs to 
include both the background noise and the photon noise 

Ktalf={<f+Hhf- (2-5) 

However, not knowing the photon noise level, we estimate the Te by setting = cr^^. In this 
paper, we estimate the background noise of the /th channel as 


oL = 5\ N!,,At, 


( 2 . 6 ) 


where 5' is the standard deviation of the data in the blue shaded region shown in Figure |2(b)| . 
As shown in Figure |2(b)| , we first fit a Gaussian function to the Thomson signal which finds the 
amplitude, peak-time and width (w) of the Gaussian. The blue shaded region where we obtain the 
background noise level starts at 5w away from the peak-time assuming that no Thomson signal 
exists in this region, hence only the background signal. is the number of data points within the 
fitted Gaussian, i.e., = 4w/A/, where At is 0.8 ns, the time step of the data points with 1.25 

GS/s. Figure ^ shows examples of the temporal evolution of Te for two different KSTAR plasma 
shots estimated using Eq. 
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Figure 4. Examples of estimated Te using a weighted look-up table method for (a) KSTAR shot #10401 and 
(b) #10433. 


3. Synthetic Thomson data and signal-to-noise ratio 


As the estimated Te shows about 15% variation on Te with respect to its mean value (e.g. see 
Figure we raise the following questions: 1) can the background noise explain such an observed 
variation on Tel 2) if not, how large does have to be to explain the observation? To be able to 
answer these questions, we generate synthetic Thomson signal where we can vary the noise levels 
numerically. 

By letting niaser (0 Eq. (pT| ) have a Gaussian form in a time domain (consistent with fitting 
the Gaussian function to measured Thomson signals), we generate the synthetic Thomson data for 
each channel similar to the measured data shown in Figure 0. We set the time step of synthetic 
data to be 0.8 ns corresponding to 1.25 GS/s of the fast sampling digitizer. A random number 
selected from a normal distribution is, then, added to each time point as a noise. The width of this 
normal distribution 8syn is set such that Gsyn = 
level. This completes generating synthetic Thomson signal for a single laser pulse. Let us denote 
the /th channel synthetic Thomson signal as V^yn (0 ii^ integral as A^^yn estimated by fitting a 
Gaussian function to the V^yn (0 if ii were real measured signal. 

We define the signal-to-noise ratio (SNR) of SNRt^, and that of the measured (synthetic) 
Thomson signal from the /th channel, SNR^^ (SNR^^yn^^ 


(2.6)) controlling the noise 


SNRt, 





SNR^syn — 



(3.1) 


where the subscript * takes any one of the total (‘total’), background (‘bg’) or photon (‘ph’) for the 
noise source, i.e., SNR^^ indicates the SNR of the Thomson signal due to the background noise. 
Here, (•) indicates the time average with 50 consecutive data points equivalent to the time duration 
of 0.5 seconds, a couple of equilibrium evolution time scale for a typical KSTAR plasma, with 
the laser pulse repetition rate of 100 Hz during the 2014 KSTAR campaign. CTj; is the standard 
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Figure 5. Examples of SNRt^ (black circle) and for i = 2 (red asterisk), 3 (blue diamond) and 4 

(green triangle) for (a) KSTAR shot #10401 and (b) #10433. Here, {Te) is estimated using chs 2, 3, and 4. 


deviation of during the 0.5 seconds. Figure ^ shows time evolution of SNRt^ and SNR^^ for 
/ = 2, 3 and 4 for the plasma shots shown in Figure 0. 


4. Estimating photon noise level and photon counts 


We wish to find and the corresponding for / = 2,3 and 4 consistent with the ex¬ 
perimentally measured SNRj^. Once we have can be calculated using Eq. (2.5) with 


the measured ^bg- Thus, we obtain 


and gL.i = CT. 


SYN 


^total 


^SYN' 


total ’ ph _ 

and with the synthetic data by recognizing that 
Here, we face a problem: there may exist infinite number 
of solutions on the combination of SNR\^^^i for / = 2,3 and 4 consistent with the measured SNRj^. 
To circumvent such a problem, anticipated by the experimental observation, we generate synthetic 
data only for chs 2 and 3 and let SNR}^^ = SNR^^^ 


since SNR^^ < SNR^^ ^ Figure 

We generate a database of SNRj^ = / {SNR^^yn^ for many different values of (Te) for i = 2 and 
3. Figure ^ shows the level of SNRp^ estimated with the synthetic data as a function of SNR^^yn 
{T e) = 800 eV. 


As the SNRp^ from the synthetic data has been estimated with only chs 2 and 3 (e.g. Figure 
we recalculate the experimental SNRp^ using only these two channels as shown in Figure 0. 
Here, we also plot the ‘required’ level of SNR\^^^i (green square) using the database of SNRp^ — 
f (SNR^^yn'> which allows us to estimate for i — 2 and 3 using Eq. (^T|). 

Figure ^ shows, for the KSTAR #10433, the levels of total required noise obtained by the syn¬ 
thetic data, background noise from experimental data and photon noise with Eq. (^31) normalized 
to the experimental signal levels for chs 2 and 3. It is clear that the KSTAR Thomson scattering 
system during 2014 campaign is photon noise dominated at least for the chs 2 and 3 of the polychro- 
mator we have investigated. This is good that the system is not limited by the background noise, 
however it also suggests that the system needs to collect more photons and utilize the collected 
photons more efficiently to reduce the 15% variation on the measured Te. 
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Figure 6. With synthetic data of chs 2 and 3, SNRt^ is estimated as a function of SNRsyn at (Te) = 800 eV. 
This allows us to obtain the required values of SNR^^yn corresponding (y^YN the experimental 

SNRt,. 




Figure 7. Experimentally measured SNRt^ (black circle) and SNR^^ for i = 2 (red asterisk) and 3 (blue 
diamond) for (a) KSTAR shot #10401 and (b) #10433. Green squares show the ‘required’ SNR\^^^i estimated 
with the synthetic data consistent with the observed SNRt^. Here, (r^) is estimated using chs 2 and 3 only. 


From the estimated normalized photon noise level, i.e., l/SNR^^^, we can estimate the ‘effec¬ 
tive’ photon counts detected by the photon detectors as 


1 _ 1 



(4.1) 


where F is the noise factor defined as the ratio of the input SNR to the output SNR 0], and Np^ is 
the ‘actual’ photon counts. Using the values from Figure @ and Eq. (P~T|) , we find that the effective 
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Figure 8. Total (black circle), background (red asterisk) and photon (blue diamond) noise levels normalized 
to the signal levels of (a) ch 2 and (b) ch 3 for KSTAR plasma shot #10433. 


photon counts are approximately 170 for chs 2 and 3. Note that the noise factor F is always greater 
than or equal to 1, thus actual photon counts may well be larger than 170, i.e., by a factor of F^. 
However, this cannot be regarded as the gain in photon counts since F decreases the SNR. 

The core KSTAR 2014 TS system has an effective F/# of 6.7 which corresponds to 17.5 msr. 
With the incident laser wavelength of 1064 nm and energy of 2.0 J, the incident photon number 
is estimated to be 1.1 x 10^^. Having approximate values of effective quantum efficiency of 15%, 
scattering length of 10 mm, optical transmittance of 40%, filter transmittance of 70%, we estimate 
photon budget to be approximately 6,200 photons per 1.0 x 10^^ m“^ of electron density. For 
Te = 800 eV, we find that the optical band-pass filters for chs 2 and 3 covers about 30% of the full 
spectrum (see Figure [T|). Thus, the expected photon counts on chs 2 and 3 are approximately 1,800 
about an order of magnitude higher than the estimated effective photon counts. This means that the 
KSTAR TS system has large potential to be improved. 


5. Conclusion 

In this paper, we estimated the time evolution of electron temperature by using the weighted look¬ 
up table method and found a 15% variation on Tg. By using the synthetic Thomson data, we have 
found that such a variation on electron temperature cannot be explained solely by the background 
noise and thus estimated the required uncertainty level consistent with the experimental obser¬ 
vation. Our study indicates that the chs 2 and 3 of the polychromator we have investigated are 
dominated by the photon noise with the average effective photon counts per laser pulse of 170. 
This suggests that improving the photon collection system and decreasing the noise factor of the 
system is required to decrease the observed scatters in electron temperature. 
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